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Starting from a Skyrme interaction we study the properties of the low-energy spectrum of 
quadrupole excitations in 90,92 Zr, 92 ' 94 Mo. The coupling between one- and two-phonon terms in 
the wave functions of excited states are taken into account. We use the finite-rank separable ap- 
proximation which enables one to perform the QRPA calculations in very large two-quasiparticle 
spaces. Our results from the SGII Skyrme interaction in connection with the volume pairing inter- 
action are in reasonable agreement with experimental data. In particular, we present the successful 
description of the Ml transition between low-energy quadrupole excitations. 



PACS numbers: 21.60.Jz, 23.20.-g, 21.60.Ev, 27.60.+J 
I. INTRODUCTION 

Properties of the quadrupolc-collective isovector exci- 
tations of the valence shell of heavy nuclei reflect three 
main aspects: collectivity, shell structure, and the isospin 
degree of freedom. The balance between these aspects 
has been studied within a series of experiments on so- 
called mixed-symmetry (MS) states. There is a rather 
complete list of references on that subject in review pQ. 
These isovector excitations were predicted in the proton- 
neutron version of the interacting boson model (IBM- 
2) [2], where the proton-neutron symmetry of the wave 
functions is quantified by the bosonic analog of isospin, F 
spin [3H5]. In particular, there are the symmetric states 
with maximum F-spin (F = F max ) and the MS ones with 
F < F max . One of the well-known examples is the MS 
quadrupole state with F = F max — 1 in the nucleus 94 Mo 
which was indentified by the measured strong Ml tran- 
sitions and weakly-collective E2 transitions to low-lying 
symmetric states [6HH]. In contrast to the case for 94 Mo, 
the next lighter even-even N = 52 isotone, 92 Zr has the 
proton subshell closure Z — 40 and the collectivity of the 
2\ state decreases. As a result, the collective MS struc- 
ture and the structures showing the considerable F-spin 
breaking have been observed experimentally [51-11 1] . This 
sensitivity makes the MS states highly appealing objects 
for studies in microscopic approaches. 

The results of the large-scale shell-model calcula- 
tions [9] [lOl O [13] have indicated the occurrence of 
MS quadrupole states in 94 Mo, 92 Zr and the breaking 
of F-spin symmetry in 92 Zr. For the case of 94 Mo, the 
wave functions, the Ml and E2 matrix elements support 
the MS assignments for the 2% state and the consistency 
with the IBM-2 predictions proves the collective charac- 
ter of the observed low-lying levels. For 92 Zr, the cal- 
culated 2f 2 states are predominantly the pure neutron 
and isovector excitations, respectively. This means that 
the single particle and collective degrees of freedom are 
present in the low-energy spectrum. 

An alternative powerful microscopic approach is the 



quasiparticle-phonon model (QPM) |14| . The model 
Hamiltonian is diagonalized in a space spanned by states 
composed of one, two and three phonons which are 
generated in quasiparticle random phase approximation 
(QRPA). The separable form of the residual interaction 
is the practical advantage of the QPM which allows one 
to perform the structure calculations in large configura- 
tional spaces. For 94 Mo, the QPM confirms the IBM-2 
level scheme, selection rules and, in particular, the dom- 
inant one-phonon structure of the transitions to the first 
and third 2+ states [HI US]. Also, the QPM gives a sat- 
isfactory and comprehensive description of the large va- 
riety of data measured for 92 Zr [TT1 OS HZj- The QPM 
wave functions of the two lowest 2 + states are dominated 
by the lowest neutron and proton two-quasiparticle com- 
ponents and the structures are very similar to results of 
the shell model calculations. Thus, the QPM can provide 
a microscopic support to the IBM-2 scheme. However, it 
is difficult to extrapolate the model parameters to new 
regions of nuclei. 

The properties of MS states are particularly sensitive 
to the proton-neutron interaction [TS1 |Tj5] . On the other 
hand, the dominance of the proton-neutron attraction is 
one of the main characteristics of the effective nucleon- 
nucleon interaction. It can be traced back to the coop- 
eration of its T = and T = 1 channels. This is a good 
possibility to examine microscopic approaches using ef- 
fective nucleon-nucleon interactions. One of the success- 
ful tools for nuclear structure studies is the QRPA with 
the Skyrme interaction [5DJ [2T] . These QRPA calcula- 
tions allow to relate the properties of the ground states 
and excited states through the same energy density func- 
tional. Although such an approach describes the proper- 
ties of the low-lying states less accurately than more phe- 
nomenological ones, the results are still in a reasonable 
agreement with experimental data, in particular, with re- 
spect to qualitative features of the properties of the 2f 
states of predominantly one-phonon excitations, see for 
example [22ti25] , 

The complexity of calculations taking into account the 
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coupling between one-phonon and more complex states 
increases rapidly with the size of the configurational 
space. Making use of the finite rank separable approxi- 
mation (FRSA) [521 HZ] for the residual interaction en- 
ables one to perform QRPA calculations in very large 
two-quasiparticle spaces. Taking into account the basic 
QPM ideas, the approach has been generalized to take 
into account the coupling between one- and two-phonon 
components of the wave functions [21] • The FRSA has 
been used to study the electric low-lying states and giant 
resonances within the QRPA and beyond [2^H50] . Alter- 
native schemes to factorize the residual interaction have 
also been considered in Refs. [3"TH3"3"] . 

Before to investigate the occurrence of MS states in 
neutron-rich nuclei one needs to be sure that the ap- 
proach is sufficiently good to reproduce characteristics of 
the low-energy spectrum of quadrupole excitations of nu- 
clei in the mass range A « 90, in particular, the Ml tran- 
sitions between the excitations. This initial paper gives 
an illustration of our approach for 92 Zr, 94 Mo in com- 
parison to the N — 50 isotones 90 Zr, 92 Mo with closed 
neutron shell. Preliminary results of our studies for 94 Mo 
are reported already in Ref. [53] . 

This paper is organized as follows. In Sec. II, we sketch 
the method allowing to consider the phonon-phonon cou- 
pling. In particular, the QRPA equations in the case 
of the finite rank separable form of the residual interac- 
tion are discussed in Sec. IIA. The coupling between one- 
and two-phonon terms in the wave functions of excited 
states are taken into account in Sec. IIB. In Sec. Ill, 
we discuss the details of our calculations and show how 
this approach can be applied to treat the proton-neutron 
mixed-symmetry states. Results of our calculations for 
properties of the quadrupole states in 90,92 Zr, 92,94 Mo 
are given in Sec.IV. Sec. IVA is devoted to QRPA analy- 
sis, while the effects of the phonon-phonon coupling are 
discussed in Sec. IVB. Conclusions are finally drawn in 
Sec. V. 



II. THE METHOD 

This method has already been introduced in Refs. [2l)l - 
129] ; hence, let us briefly summarize the different steps. 
For the present study, the anharmonicity of the low- 
energy vibrations is constrained by the coupling between 
one- and two-phonon terms in the wave functions of ex- 
cited states. 



A. Implementation of QRPA 

The starting point of the method is the HF-BCS calcu- 
lation 35 of the ground state, where spherical symmetry 
is imposed on the quasiparticle wave functions. The con- 
tinuous part of the single-particle spectrum is discretized 
by diagonalizing the Skyrme HF Hamiltonian on a har- 
monic oscillator basis. As effective interactions, Skyrme 



interaction are used in the particle-hole (p-h) channel, 
and the pairing correlations are generated by a density- 
dependent zero-range force 

V p ^r{Ti,Ta) = V (l-ri(^^j )*(ri-r 2 ), (1) 

where p{r\) is the particle density in coordinate space; 
Po is equal to the nuclear saturation density; a, r\ and Vq 
are model parameters. We use a=l; r\— 0, r\— 0.5, and 
rj= 1 are the cases of a volume interaction, a mixed in- 
teraction and a surface-peaked interaction, respectively. 
The strength Vq is fitted to reproduce the odd-even mass 
difference in the region of nuclei considered here. In or- 
der to limit the pairing single-particle space, we have 
used the smooth cutoff at 10 MeV above the Fermi ener- 
gies EH EH Eg . 

The residual interaction in the p-h channel and 
in the p-p channel V pp s can be obtained as the second 
derivative of the energy density functional with respect 
to the particle density and the pair density, accordingly. 
As proposed in Ref. [26] , we simplify V p ^ s by approximat- 
ing it by its Landau-Migdal form. All Landau parame- 
ters with I > 1 are equal to zero in the case of Skyrme 
interactions. We keep only the I = terms in V p £~ s - In 
this work we study only normal parity states and one 
can neglect the spin-spin terms since they play a minor 
role [35] [25] ■ Also, the two-body Coulomb and spin-orbit 
residual interactions are dropped. Therefore we can write 
the residual interaction as 

V; a es ( ri ,r 2 ) = N - 1 [FS(r 1 )+F^(r 1 ) Tl -T 2 ]d(ri~r 2 ) (2) 

where a is the channel index a — {ph,pp}; t$ are the 
isospin operators, and No = Ikpm* /ir 2 h 2 with kp and 
m* standing for the Fermi momentum and nucleon effec- 
tive mass. 

The expressions for FQ h ,F Q ph and F PP ,F Q PP can be 
found in Ref. [35] and in Ref. [571, respectively. Since 
the definition of the pairing force (pj involves the energy 
cutoff of the single-particle space to restrict the active 
pairing space within the mean-field approximation, this 
cutoff is still required to eliminate the p-p matrix ele- 
ments of the residual interaction in the case of the sub- 
shells that are far from the Fermi energies as in Ref. [37] . 

The p-h matrix elements and the antisymmetrized p-p 
matrix elements can be written as the separable form in 
the angular coordinates (26] [27] . After integrating over 
the angular variables we use a iV-point integration Gauss 
formula for the radial integrals. Thus, the residual inter- 
action can be expressed as the sum of N terms in FRSA 
for the Skyrme residual interaction (26] [27] . 

We introduce the phonon creation operators 

jj' 

-(-l) A -^U(jj';A-/i)), (3) 
A + (jf; A/i) = ^ (jmj'm' \ Xfi)aj m a+ m ,, (4) 

mm/ 
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where A denotes the total angular momentum and /i is 
its z-projection in the laboratory system. a>j m (ctj m ) is 
the quasiparticle creation (annihilation) operator and jm 
denote the quantum numbers nljm. One assumes that 
the ground state is the QRPA phonon vacuum | 0) and 
the one-phonon excited states are | 0) with the nor- 
malization condition 
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Making use of the linearized equation-of-motion ap- 
proach one can get the QRPA equations [35]. Solutions 
of this set of linear equations yield the one-phonon en- 
ergies u>\i and the amplitudes Xfy, Y^, of the excited 
states. There are two types of QRPA matrix elements: 



the A 
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matrix related to forward-going graphs 



and the B 
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'a ;i\u matrix related to backward-going 
graphs, and the dimension of these matrices is the space 
size of the two-quasiparticle configurations. Using the 
FRSA for the residual interaction, the eigenvalues of the 
QRPA equations can be obtained as the roots of a rela- 
tively simple secular equation, where the matrix dimen- 
sions never exceed 6iV x 6iV independently of the two- 
quasiparticle configuration space size [5S1 [37] . The stud- 
ies [2"TH2"9"] enable us to conclude that 7V=45 for the rank 
of our separable approximation is enough for the electric 
excitations (J < 6) in nuclei with A < 208. 



B. Phonon-phonon coupling 

In the next stage, we construct the wave functions from 
a linear combination of one-phonon and two-phonon con- 
figurations 



AlilA2»2 

with the normalization condition 
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Because of the completeness and orthogonality condi- 
tions for the phonon operators, the bifermion opera- 
tors A + (jj'; Xfi) and A(jj'; A/z) can be expressed by the 
phonon ones. The Hamiltonian can be rewritten in terms 
of quasiparticle and phonon operators [141 128) . 

The amplitudes Ri{Xv) and P^(Xv) arc determined 
from the variational principle 



«*„(Aa*)|*„(A/x)> - 1)) = 0, 



(8) 



which leads to a set of linear equations [T3] |2"5] 
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■U^(\ l )R l (Xiy) = 0. 



(10) 



The Pauli principle corrections are dropped here since 
the effect on the lowest excited states is small 14J. 

The rank of the set of linear equations ([9]), ( 10 1 is equal 
to the number of one- and two-phonon configurations in- 
cluded in the wave function (|6j). Its solution requires 
to compute the matrix elements coupling one- and two- 
phonon configurations 



(0\QmH [Qt lh Q 
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Eqs. (9|, (10) have the same form as the QPM equa- 



tions [14], but the single-particle spectrum and the pa- 
rameters of the residual interaction are calculated with 
the Skyrme forces. 



III. DETAILS OF CALCULATIONS 

We apply the approach to study the low-energy spec- 
trum of 2 + excitations in the nuclei 90,92 Zr, 92,94 Mo. 
In our calculations the single-particle continuum is dis- 
cretized by diagonalizing the HF Hamiltonian on a basis 
of 12 harmonic oscillator shells and cutting off the single- 
particle spectra at the energy of 100 MeV. This is suffi- 
cient to exhaust practically all the energy-weighted sum 
rule. Because of the large configurational space, we do 
not use effective charges. 

The single-particle structure around the Fermi level 
is a key ingredient in the microscopic analysis. For 
the investigation, we adopt five Skyrme interactions, 
namely, SGII [35], SLy4 [31], SLy5 [3J], SLy5+T 00], 
and T44 [IT]. The SGII parameterization is the suc- 
cessful attempt to describe the spin-dependent proper- 
ties from a standard Skyrme force. In particular, one ob- 
tains a good description of experimental energies of the 
Gamow-Teller resonance of 90 Zr [35]. The parameters of 
the forces SLy4 and SLy5 have been adjusted to repro- 
duce nuclear matter properties, as well as nuclear charge 
radii, binding energies of doubly magic nuclei. The forces 
SLy5+T and T44 involve the tensor terms added with 
(T44) and without (SLy5+T) refitting the parameters of 
the central interaction. These five parametrizations de- 
scribe correctly the subshell order near the Fermi level of 
90 ' 92 Zr, 92 ' 94 Mo. To see this, the calculated neutron and 
proton single-particle energies for the case of 92 Zr and the 
experimental data [H] are plotted in Fig. 1. The calcu- 
lations with SGII reproduce well the experimental data 
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FIG. 1: Neutron (top) and proton (bottom) single-particle energies (in MeV) near the Fermi energies for 92 Zr calculated with 
SGII, SLy4, SLy5, SLy5+T, and T44. The experimental spectra are taken from Ref. [42]. 
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FIG. 2: Energies (in MeV) of lowest two-quasiparticle states in 92 Zr. The results of calculations with the volume, the mixed, 
and the surface-peaked pairing interactions are shown. 



and other choices of the Skyrme forces do not improve 
the agreement. 



We employ the isospin-invariant pairing force ([I]) , with 
the value /?o=0.16 fm -3 of the nuclear saturation density 
corresponding to the SGII force. The strength Vq is fit- 
ted to get a reasonable description of the experimental 



neutron and proton pairing energies of 90,92 Zr, 92,94 Mo, 
Pn = \{B (N, Z) + B (N - %Z) - 2B {N - 1, Zj) 

for neutrons, and similarly for protons. Thus, the 
strength V is taken equal to -270 MeVfm 3 , -420 MeVfm 3 , 
and -870 MeVfm 3 for the cases of the volume, the mixed 
and the surface-peaked interaction, respectively. In or- 
der to make the choice of the pairing interaction, it is 
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FIG. 3: Energies and B(E2)-values of the lowest quadrupole 
states in 92 Zr. The QRPA calculations are performed without 
(dashed line) and with (solid line) the residual p-p interaction. 



very useful to analyze the three lowest two-quasiparticle 
states in 92 Zr. Since the 2d^/ 2 neutron subshell is par- 
tially filled and the 2pi/ 2 proton subshell is filled, one can 
expect that the first level is the {2^5/2,2^5/2} neutron 
state while the second level is the {139/2,139/2} proton 
one, see for example [TJ [TBI [T7] . As it is seen from Fig. 2 
the different choices have a strong impact on the unper- 
turbed two-quasiparticle excitations. Only the volume 
pairing interaction gives the {159/2,139/2} proton state 
as the second level. Thus, hereafter we use the Skyrme in- 
teraction SGI1 in the particle-hole channel together with 
the volume zero-range force acting in the particle-particle 
channel. 

The Landau parameters Fq H , Fg h expressed in terms 
of the Skyrme force parameters [35] depend on ^. As 
it is pointed out in our previous works US] one needs 
to adopt some effective value for kp to give an accurate 
representation of the original p-h Skyrme interaction. For 
the present calculations we use the value k F — I.45fm _1 
adjusted so as to have the spurious center-of-mass state 
at zero energy. 

It is worth to mention the effect of the residual p-p 
interaction within the QRPA. Fig. 3 shows the effect on 
the 2^ energies and the B(E2) values of 92 Zr. Including 
the quadrupole p-p interaction results in a decrease of the 
2^ energies (about 0.1 MeV) and the -B(£?2)-values do 
not change practically. 

Finally, we discuss the extension of the space for one- 
and two-phonon configurations. To construct the wave 
functions (6j of the low- lying 2 + states up to 4 MeV we 
use only the 2+ phonons, and all one- and two-phonon 
configurations with energies up to 8 MeV are included. 
We have checked that the inclusion of the high-energy 
configurations plays a minor role in our calculations. 



IV. RESULTS 



A. QRPA analysis 



First, properties of the low-lying quadrupole states are 
studied within the one-phonon approximation. Results 
of our calculations for the 2+ 2 states in 90 ' 92 Zr, 92,94 Mo: 
the energies, the B(E2), B(M1) values, and the con- 
tributions greater than 10% to the normalization condi- 
tions |5| are given in Table I. Note that the B(M1) values 
have been calculated with the 3-factors of free protons 
and neutrons. 

One can see that we find a satisfactory description of 
the isotopic dependence of the 2+ energies near closed 
shells. The closure of the neutron subshell I39/2 in 90 Zr, 
92 Mo leads to the vanishing of the neutron pairing and 
as a result energies of the first two-quasiparticle poles, 
{lff9/2, 139/2}^ in 90 Zr, 92 Mo are larger than energies 
of the first poles, {2d 5/2 , 2d 5/2 }„ in 92 Zr, 94 Mo. This 
yields that the 2f state has collective (noncollective) 
structure with the domination of the proton configura- 
tion {I39/2, 139/2} for the case of 90 Zr ( 92 Mo). On the 
other hand, in 92 Zr the leading neutron configuration 
{2d 5 / 2 , 2d 5 / 2 } gives a contribution of 79% that is almost 
twice larger than in 94 Mo. The structure peculiarities 
are reflected in the B(E2: 0+ — > 2^ ) values, as is shown 
in Table I. The 2\ states in 92 Zr, 94 Mo have an isoscalar 
character since the dominant neutron and proton phonon 
amplitudes X, Y of the 2^ state are in phase. 

There are the second 2 + states in 92 Zr, 94 Mo below 
3 MeV within the one-phonon approximation. The dom- 
inant neutron and proton amplitudes of the fairly collec- 
tive 2 2 state are out of phase. As a consequence, the 
isovector character of the 2 2 states are reflected in the 
remarkable values of B(M1; 2 2 — )• 2+), as given in Table 
I. The first calculations of the isovector behavior of 2 2 
QRPA excitations in 92 Zr, 94 Mo based on the analysis of 
the phonon amplitudes within the QPM have been done 
in |I5HIZ|. 

We turn now to the structures of the 2 2 states of 
90 Zr, 92 Mo. The proton configurations exhaust about 
99% and 56% of the wave function normalization in 90 Zr 
and 92 Mo, respectively. It means that the second pole, 
{2pi/ 2 , 2p 3 f 2 } 7r in 92 Mo is closer to the neutron poles 
than that in 90 Zr. As expected, the negligible size of 
the B(M1;2+ -> 2+) value of 90 Zr is obtained. The 2+ 
state in 92 Mo has the main neutron and proton phonon 
amplitudes in phase (Table I) and this results in the com- 
parable B(M1) value of the Ml transitions between the 
noncollective 2f state and the isoscalar 2 2 state. 

This analysis within the one-phonon approximation 
can help to identify the MS state, but it is only a rough es- 
timate. Some overestimate of the experimental energies 
(Table II) indicates that there is room for two-phonon 
effects. 
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TABLE I: Energies, transition probabilities, and structures of the QRPA quadrupole states in 90 - 92 Zr and 92,94 Mo. Phonon 
amplitude contributions greater than 10% are given. 





State 


Energy 
(MeV) 


B{M1;2+ -> 2+) 


B(E2;0+ -+2+) 
(e 2 fm 4 ) 


{mhji, mhj2} T 


X 


Y 


% 


90 Zr 




2.8 




630 


{2d 5 / 2 , 109/2}* 

{109/2, 109/ 2 }tt 

{2pi/ 2 , 2p 3 /2K 


-0.37 
1.03 
-0.52 


-0.11 
0.06 
-0.03 


13 
53 
26 




2 2 + 


3.4 


0.00 


10 


{2p 1/2 ,2p 3/2 }^ 

{109/2, l09/ 2 }?r 


0.79 
0.85 


0.00 
0.00 


63 
36 


92 Zr 




1.7 




410 


{2d 5 / 2 , 2d b / 2 }v 


1.26 


0.12 


79 






2.9 


0.53 


310 


{2d 5 / 2 , 2d 5 / 2 }„ 
{3si/2, 2d 5 / 2 }i/ 
{109/2, I09/2K 

{2pi /2 ,2p 3/ 2}^ 


-0.63 
-0.45 
0.85 
-0.36 


0.09 
-0.04 
0.04 
-0.02 


20 
20 
36 
13 


92 Mo 


n 


1.9 




1170 


{2d 5/2 , 139/2};, 
{109/2, I09/2K 


-0.35 
1.32 


-0.16 
0.19 


10 

86 




2 2 + 


4.4 


0.25 


230 


{2d 5/2 , l0 9 / 2 }„ 
{2p 1/2 ,2p3 /2 }^ 

{109/2, 109/2/tt 


-0.65 
-0.63 
-0.49 


-0.07 
-0.02 
0.10 


42 
40 
11 


94 Mo 




1.2 




1730 


{2d 5 / 2 , 2d s / 2 }„ 

{109/2, 109/2K 


0.92 
0.99 


0.27 
0.37 


39 
42 




2t 


2.4 


1.23 


160 


{2d 5 / 2 , 2d5/ 2 }i/ 

{109/2, 109/2}^ 


-1.08 
0.89 


0.07 
0.02 


58 
40 



TABLE II: Energies, transition probabilities and dominant components of phonon structures of the low-lying quadrupole states 
in 90 - 92 Zr and 92 94 Mo. Experimental data are taken from Refs. P [71 ITOl l4"3l [44] . 





Af = 2+ 


Energy 
(MeV) 
Expt. Theory 




Structure 


B(E2;Q+ ->2+) 

(e 2 fm 4 ) 
Expt. Theory 


B(M1;2+ - 
Expt. 


,2+) 
Theory 


90 Zr 


2t 


2.186 


2.6 




93%[2+] QRPA 


643±22 


600 








2 2 + 


3.308 


3.2 




95%[2+] QRPA 


53±14 


1 


0.088±0.025 


0.00 


92 Zr 




0.934 


1.6 




96%[2+] QRPA 


790±62 


420 








2+ 

Z 2 


1.847 


2.7 




87%[2+] QRPA 


419±49 


230 


0.37±0.04 


0.41 




n 


2.067 


2.6 


45% [2+ 


]qapa+37%[2+ ® 2+]q RPA 


< 0.62 


50 


<0.024 


0.17 


92 Mo 




1.509 


1.9 




99%[2+] QRPA 


1036±62 


1160 








n 


3.091 


3.8 




91%[2+ <g> 2+] QRPA 


254±20 


50 


0.043±0.007 


0.03 


94 Mo 




0.871 


0.5 




73%[2+} QRPA 


2031±25 


1280 








2+ 

Z 2 


1.864 


1.8 


53% [2+ 


® 2+] QR pa+21%[2+] QR pa 


32±7 


170 


0.06±0.02 


0.07 




2+ 


2.067 


2.3 




87%[2t} QR p A 


279±25 


310 


0.56±0.05 


0.68 
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B. Effects of the phonon-phonon coupling 

Let us now discuss the extension of the space to one- 
and two-phonon configurations. The calculated 2 + state 
energies, the largest contributions of the wave function 
normalization (frfc, the B(E2) 7 B(M1) values and the ex- 
perimental data0[2[lOl|43lll4] are shown in Table II. It 
is worth pointing out that we get the wrong energy order 
of the second and third 2 + states in the case of 92 Zr. On 
the other hand, the energy difference equal to 0.1 MeV 
is close to the expected accuracy of our calculations. 

One can see that the inclusion of the two-phonon terms 
results in a decrease of the 2^ energies and in a reduc- 
tion of the B(E2; 0+ -> 2+) values, except for 92 Zr. Our 
calculations reproduce well a general behavior for ener- 
gies and transition probabilities. There is some under- 
estimation of the B(E2) values of the N=52 isotones in 
comparison with the experimental data, this probably 
points to a particular problem due to the effective in- 
teraction rather than to a deficiency of our variational 
space. In all four nuclei, the crucial contribution in the 
wave function structure of the first 2 + state comes from 
the [2f ]qrpa configuration, but the two-phonon contri- 
butions are appreciable. This means that the structures 
of the first 2 + states do not change practically due to the 
effects of the phonon-phonon coupling. As a result, we 
get the isoscalar collective structure of the 2f state in 
94 Mo and the neutron dominated 2\ excitation of 92 Zr 
which indicates the F-spin breaking. 

The second 2 + state contains a dominant two-phonon 
configuration [2+ ®2 j [\qb.pa in 92,94 Mo and such contri- 
bution leads to the small B(E2) values. Moreover, the 
calculated B(M 1) value is sensitive to the phonon com- 
position: the two-phonon configuration composed of the 
isoscalar collective phonons in the case 94 Mo and the non- 
collective phonons of 92 Mo. In 90,92 Zr, the wave function 
of the 2 2 state is dominated by the [2t\QRPA configu- 
ration and we can follow the QRPA estimate discussed 
in Sec. IVA. For the case of 92 Zr we obtain the isovector 
collective 2^ state and the dominant one-phonon struc- 
ture of the Ml transition between the first and second 
2 + states in 92 Zr. The results are close to those that were 
previously calculated in the shell model H3] and in 

the qpm nana mi- 

Finally, we examine the occurrence at low energy (be- 
low 3 MeV) of the third 2 + states. One can see that 
the collective state in 94 Mo is dominated by the isovector 
one-phonon Q^QRPA structure. The calculated values of 
the 2^ energy and the transition probabilities are in rea- 
sonable agreement with the experimental data. In other 
words, we reproduce the IBM-2 level scheme and the cal- 
culated B(M1; 2^ —> 2+ ) value supports the MS assign- 
ments observed experimentally and theoretically for the 
first time in Ref. [BJ. It is noteworthy that this conclu- 
sion for 94 Mo remains valid for the SLy5+T parameter 
set [34]. For the 2^ state of 92 Zr, one of the main com- 
ponents of the wave function is the [2j~ ® 2f]QppA con- 
figuration. As can be seen from Table II, the two-phonon 



contribution is reflected in the small values of the tran- 
sition probabilities, but there is some overestimation in 
comparison with the experimental data. One can expect 
an improvement if the variational space is enlarged by 
the phonon composition with hexadecapole multipolar- 
ity and the three-phonon configurations are taken into 
account. Such calculations are now in progress. 

Thus the microscopic approach [26H28] describes the 
properties of the low-lying states in 92 Zr, 94 Mo less accu- 
rately than more phenomenological ones [5J H"5Tfl"7] , but 
the results are still in a reasonable agreement with the 
experimental data [6l 171 HOl 143] 144] . 



V. CONCLUSIONS 

Starting from the Skyrme mean-field calculations, we 
have studied the effects of the phonon-phonon coupling 
on the properties of the low-energy spectrum of 2 + exci- 
tations and, in particular, on the Ml transitions between 
the excited states of nuclei in the mass range A ss 90. The 
finite-rank separable approach for the QRPA enables one 
to perform the calculations in very large configurational 
spaces. 

The parametrization SGII of the Skyrme interaction is 
used for all calculations in connection with the volume 
zero-range pairing interaction. Using the same set of pa- 
rameters we have studied the behaviour of the energies, 
the B(E2; 0+ 2+) and B(M1; 2+ 2+) values of the 
lowest 2 + states in 90,92 Zr, 92,94 Mo. Among our initial 
motivation was the search for MS states in 92 Zr, Mo 
in comparison to the N = 50 isotones 90 Zr, 92 Mo with 
closed neutron shell. Our results indicate indeed the oc- 
currence of MS states in our calculation for the nuclei 
92 Zr and 94 Mo that were successfully predicted within 
the IBM-2 before. Our results from the Skyrme interac- 
tion are in reasonable agreement with experimental data. 
We stress that they represent the first successful compari- 
son between experimental Ml transition values and those 
calculated with the Skyrme interaction. The coupling be- 
tween one- and two-phonon terms in the wave functions 
of excited states is essential. The QRPA results tend 
to overestimate the 2f energies and the inclusion of the 
two-phonon configurations results in a decrease of the 
energies. There is a clear influence on the structure of 
the 2J3 states. The structure of the low-lying 2 + states 
calculated in our approach arc close to those that were 
calculated within the QPM before. We conclude that the 
present approach may provide a valuable globally appli- 
cable microscopic analysis of the properties of the lowest 
quadrupole excitations. 

Our model would probably be improved by enlarging 
the variational space for the 2 + states with the inclu- 
sion of the two-phonon configurations constructed from 
phonons with hexadecapole multipolarity and taking into 
account the three-phonon configurations. The computa- 
tional developments that would allow us to conclude on 
this point are still underway. 
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